SESSION I: THE NATURE OF THE DATA 


Julius London presented a "tutorial overview" which began with a display 
the locations of total ozone stations (Fig. 1) and of stratospheric ozone samp 
lings (Fig. 2). The samplings are largely concentrated in three areas--Japan, 
Europe, and India. Approximately 75% of the total ozone measurements are made 
with Dobson instruments which offer the best international measurements. If th 
are well cared for and well calibrated, their accuracy is on the order of a 
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Figure 1. Location of Total Ozone Stations (1957-1975) 



Figure 2. Observations of Stratospheric Ozone (1957-1975) 

(non-satel 1 ite) 


8 


■ I ■■ !■ ■■■ 


II II I iiin I 


III 


II Hill llllll I lllll■■■■llllll 






few percent. Long period observations made with the same instrument have "rela- 
tively good" reproduction and reliability. The remaining 25% of the world's 
observations are made with filter instruments which were very noisy until the 
USSR improved them substantially in 1969. "Now," London said, "they are noisy 
rather than notoriously noisy." 

Plotting the available data produces a picture of global dis- 
tribution of total ozone (Fig. 3) that shows an equatorial minimum 
and an increase toward the polar regions. These features of the 
ozone distribution are both well-known. Variations in total ozone 
amount depend on season and latitude. This can also be seen in 
variance data from individual stations (Fig. 4). 



Figure 3. Distribution of Total Ozone for the Period 1957-1975 

(m atm-cm) 


Although the total ozone pattern is similar in both hemispheres, 
the northern hemisphere has 3 to 10% more total ozone than the 
southern hemisphere, not a negligible amount, according to London. 

In the stratosphere, the hemispheric difference is even more 
pronounced as mid-latitude eddy transport is stronger in the 
northern hemisphere than in the southern. The maximum ozone con- 
centration that occurs in the lower stratosphere varies with 
latitude and season. (In the summer and at the equator, the 
maximum is higher than in the winter and at the poles.) 

London noted the "close association" between total ozone dis- 
tribution and pressure distribution in the atmosphere. High amounts 
of ozone are clearly associated with large scale troughs in pres- 
sure distribution. 
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Displaying the average latitude seasonal variation (Fig. 5) 
shows the well-known spring maximum that is stronger in the 
northern hemisphere than in the southern. This spring maximum 
is also delayed a month in the south, a fact that London said was 
“more evidence of the relationship between ozone concentration 
and circulation" because it reflected the longer relaxation time 
in terms of the atmospheric circulation. 



Figure 5. Average Latitude Seasonal Variation of Total Ozone 



Figure 6. Latitudinal Distribution of the Harmonic Parameters 

of the Ozone Variation 

To further reinforce his point that the variance of total 
ozone concentration is very like the behavior of the lower 
stratosphere and therefore the result of meteorological parameters, 
London fitted the mean monthly latitude values of ozone to the 
first harmonic (annual cycle) and second harmonic (phase) of the 
annual ozone variance (Fig. 6). Four different comparisons of 
these data (Fig. 7) then further supported his thesis. 

Several statistical scientists commenting on their own efforts 
to smooth out the seasonal variance noted that: the variance of 

the reciprocals is less than that of the Dobson data; and nothing 
really stabilizes the variance although the reciprocal is better 
than the logarithm. 
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Figure 7. Comparison between the seasonal variation (monthly mean 

values) ( ), the day-to-day changes within 1 month (February 

1973), given by daily mean values ( ^), possible variation 

within 1 day (16 April 1962, single readings) ( ■), and the 

differences between mean values of the same calendar month (February) 
from year to year) ( ). 



Figure 8. Ten Year Running Trend, Arosa, 1932-1973 (monthly means) 


London then presented data from a single station, Arosa, 
Switzerland, where total ozone measurements have been made since 
1926. These data (Fig. 8) indicate a long-term change of 6%, 
about half the standard deviation. 

Walter Komhyr presented detailed information about Dobson 
measurements and the various algorithms used to calculate total 
ozone. The Dobson spectrophotometer (Fig. 9) takes in light 
from the direct sun or the zenith sky. The measurement prin- 
ciple requires measuring the ratio of intensities of solar 
radiation I/I' at pairs of wavelengths, A and A' (Fig. 10). 
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Figure 9. The Optical System 
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Figure 10. Ratio of Intensities of Solar Radiation I/I' at Pairs 

of Wavelengths, X and X' 

Total ozone x is then calculated from the following equation: 

Lq - L - (6 - 3')mp/Po - (6 - 6') sec Z 
^ ~ (a - a' )y ^ 

where L = log(I/r ); 

Lq is the value of L at the top of the atmosphere; 
a and a' are ozone absorption coefficients; 

3 and 3' are molecular scattering coefficients; 

6 and 6' are particle scattering coefficients; 
p is the slant path of light through the ozone layer; 
m is the air mass; 

p and Pq are mean station pressure and mean sea level 
pressure, respectively; 

and Z is the solar zenith angle. 

The ozone absorption coefficients used before 1956 were 
36% different from those currently set by the International 
Ozone Commission (Fig. 11). 
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Wavelength 
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Mean 
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A.U. 

Equivalent 
Instrument 
Slit Widths 
A.U 

Atmos pheri c 

Scattering Coefficients 
[Rayleigh Cabannes fbrmula) 

Ozone Absorption Coefficients . 

Adopted for ""^..Adopted for Use 
Use July 1 , 1957*^\January 1 , 1968** 

X 

x' 

X 

x' 

X 

x' 

0 

® B — b' 

a 

a 

a — a' 

a — a 

A 

Short 

Lona 

3055 

3254 

9 

30 

0.491 

n.375 

0.116 

1.882 

0.120 

1 .762 

1.748 

B 

Short 

Long 

3088 

3291 

9 

30 

0.470 

0.357 

0.113 

1.287 

0.064 

1.223 

1.140 

C 

Short 

. Lonq 

3114.5 

3324 

9 

30' 

0.453 

0.343 

0.110 

0.912 

0.047 

0.865 

0.800 

D 

Short 

Lonq 

3176 

3398 

9 

30 

0.417 

0.312 

0.104 

0.391 

0.017 

0.374 

0.360 

c' 

Short 

Long 

3324 

4536 

30 

15 

0.343 


0.047 

nil 

0.047 


■AD 





0.012 


1.388 

1.388 

BD 





0.009 


0.849 

0.780 

AC 

CD 





0.006 


0.897 

0.948 





0.006 


0.491 

0.440 


♦From July 1, 1957, the values of a are based on 1953 results of Vigroux for -44“C, 
which are about 36X smaller than the values of Ny and Choong used previously. 

♦♦Based on coefficients remeasured by Vigroux In 1967 as well as on atmospheric 
measurements, and recoimended for use by the lAflAP. 


Figure 11. Dobson Spectrophotometer Wavelengths and Constants 


Solution of equation (1) involves dealing with the usually unknown 
term (6 - 6'). In practice, therefore, observations are made on 
double pair wavelengths in which case the particle scattering 
term becomes (6 - 6 ' - (6 - 6 ' ) 2 - To a good approximation, 
this quantity is assumed to equal zero. The equation for 
computing ozone using two wavelength pairs is 


Xi2 


[(L^-L)j-(Lp-L)^-[(6-3')i-(3-3')2]mp/pQ-[(6-6')j-(5-6')2]sec Z 
“ — — [{a-ot' )j-(a-a' )2]y 


( 2 ) 


The most commonly used wavelength pairs are the A and D pairs 
and this equation is 


do - L)i - do - L )2 0.012mp 

^AD = 1 . 388y “ ypo 


(3) 
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Values for L, 6 , g', a, a', y, m, p and p are obtained 
using laboratory or other techniques. Calibration, i.e., measure- 
ments on direct sun, are needed to obtain Lq. The following 
equations describe the procedure for calibrating Dobson ozone 
spectrophotometers on an absolute scale. 


Method 1. (Langley Plots). For observations made on clear days 
when (6 - 6 ') - 0 , equation ( 1 ) can be rewritten 

L + (B - B')mp/pQ = -(a - a')yx + Lq (4) 

which is linear in y of the form y = ay + b provided that x re- 
mains constant during the observing interval. By plotting 
L + (B - B')mp/Po against y (for 0 s u S 3.2) and fitting a line 
to the data, the slope is a = -(a - a')x and the intercept is 
b = Lq. In this manner the extra-terrestrial constant, Lq, is 
determined for the instrument. 

Method 2. (Slope Method). Equation (1) can also be written 

Lq - L - (B - B')mp/pQ = (a - a')yx + (6 - 6 ') sec Z (5) 


Let L^ be the assumed approximate value of Lq, Lq be the true 
value, and define S = Lq - L$. Then the equation 


L* - L - (B - 3')mp/pQ 

y 


S sec Z 

- + (a - a')x + (6 - 6 ') ( 5 ) 

y y 


is linear in 1/y. Plotting the left hand side of equation ( 6 ) against 
1/y gives the slope a = -S, and the intercept is b = (a-a')x + (6-6') 
for sec Z - y. The estimate of S is then used to correct the 
estimated Lq value. 

The primary standard Dobson spectrophotometer (instrument #83) 
was used as a reference instrument for a World Meteorological 
Organization-sponsored International Comparison of Dobson Ozone 
Spectrophotometers Meeting held in Boulder, Colorado, August 8-19, 
1977. Regional secondary standard instruments from Australia, Canada, 
East Germany, Egypt, Japan, the United Kingdom, and India were 
compared with reference instrument #83. 


A summary of similar instrument comparisons held in the past 
was presented by Komhyr (Fig. 12). 
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1 . 

Hungary, 1969: 







5 Eastern European Dobson instruments 

compared 



2 groups insts. 

AO3 ~ 10% 





Discrepancies exceeding 

20% 




2. 

Belsk, Poland 1974*: 






Inst 

No 

— 

^^AD’ 

X =0.300 

cm 


% error in 

X 

^”^AD 

y=l 

y=2 

y=3 

y=l 

y=2 

y=3 

41 

0.016 

0.012 

0.006 

0.004 

3.8 

1.9 

1.3 

64 

0.108 

0.079 

0.039 

0.026 

26.3 

13.0 

8.7 

77 

-0.023 

-0.017 ■ 

■0.008 

0.006 

-5.5 

-2.8 

-1.8 

00 

0 

0 

0 

0 

0 

0 

0 

84 

0.056 

0.040 

0.020 

0.013 

13.5 

6.7 

4.5 

96 

0.002 

0.001 

0.001 

0 

0.5 

0.2 

0.1 

101 

0.018 

0.013 

0.006 

0.004 

4.3 

2.2 

1.4 

108 

0.025 

0.018 

0.009 

0.006 

6.0 

3.0 

2.0 

no 

0.051 

0.037 

0.018 

0.012 

12.3 

6.1 

4.1 

112 

-0.024 

-0.017 - 

■0.009 - 

•0.006 

-5.8 

-2.9 

-1.9 

*Spectrophotometer No. 83 
comparisons . 

was the 

reference instrument for 

the 


Figure 12. Results of Past Dobson Instrument Comparisons 


Komhyr then presented data on additional absolute calibrations 
of U.S. standard spectrophotometer #83 that were made in 1962 and 
1972. He also presented an analysis of the data which illustrated 
a method whereby possible variations in L values for the A, C, 
and D wavelength pairs may be detected. °Such variations, for 
example, may occur with variations in sunspot numbers. 

Komhyr concluded his presentation by examining the effect of 
possible changes in Lq with time on the accuracy of ozone measure- 
ments. He examined three cases: 
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Case 1 : Single pair wavelengths A, B, C, D 


Lq - L - (e-3') mp/po “ (5-6') sec Z 



* ' *-0 increase during the 1960's, measured 0^ 

amounts were too low. 

Assume: 

AL = 0.015 and x = 0.300 cm. 

° 1976-1962 


Then using: 
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Assuming that for the time interval 1962 to 1976 


= 0.0154 aLqp = 0.0141 x = 0.300 

. _ 0.0013 

1.388 M 


% Error in x^^ 

P = 1 

0.3 

P = 2 

0.2 

P = 3 

0.1 


Komhyr pointed out that if it is assumed that solar irradi- 
ance variations occur exponentially in the uv region of the 
solar spectrum, then errors in total ozone measurements re- 
sulting from Lq variations are insignificant when observations 
are made on double pair wavelength such as the AD pair which is 
the usual practice. However, observations made on single pair 
wavelengths could be significantly in error. 


Case 3: CC Observations 

These are ozone observations made on 
using C wavelengths. 

For ALqC 0.015 and x = 0.300 

the zenith 

% Error in x^^i 

y = 1 

4.0 

y = 2 

2.7 

y = 3 

2.0 


Thus the errors in CC observations due to solar 
intensity variations in the uv region of the solar 
spectrum can be appreciable. 



Reid Basher then discussed systematic errors in Dobson measure 
ments. He pointed out that the measurements are always in error 
to some extent, as a result of such things as aerosol scattering 
character, other atmospheric absorbing species, stratospheric 
temperature, instrument temperature, solar zenith angle, and solar 
spectral irradiance, but he noted that only those with long-term 
periodicities are important. 

Instrument-^related errors can be controlled by regular, 
rigorous calibration, but, in practice, this may be difficult, 
and not all errors will be known or appreciated. Atmosphere- 
related errors - attenuation by aerosols, absorption 
by other species such as SO2 and N02» dependence of ozone ab- 
sorption on temperature, and variations in the effective mean 
height of the ozone layer - can conceivably have long-term peri- 
odicities. For example, atmospheric particulate loading rises 
sharply with volcanic activity and decays with a relatively 
long half-life. These variable error sources will probably need 
to be monitored simultaneously in order to attain a high ozone 
trend detectability. 

Basher selected as an example the case of atmospheric aerosols 
particulates, and absorbers, explaining as follows: 

Since the measurement between a pair of wavelength bands is 
a differential one, only relative effects are important, e.g.: 


m (3^-62) - sec Z ((Sj-62) 

t ' , 

aerosol 

mol ecul es 

(Ti-Yj) 

t 

Other 


log 

t 


= log 


'ol 


T T 

^2 ^o2 

t 


- y X (a^-a2) - 

t 


measurement ozone 

extraterrestrial 
constant 


Consider the effect of aerosol attenuation. If the effect 
is spectrally flat, 6, = 62 “ < 5 (Xo)» it is eliminated directly 
within a band-pair measurement. If it is spectrally linear, 

^1 ~ ^^^0^ ^ ^l^^o ~ ^1^’ ® ^hird wavelength band or second band- 

pair is needed to eliminate the error. 
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2 

For a quadratic dependence, 6^ = <s(x^) + gj^(x^-x^) + 92 (^ 0 ”^!^ ’ 

four wavelengths or three band-pairs are needed to eliminate the 
error. The standard Dobson measurement, X«q, is a two band-pair 
measurement and thus deals effectively witn spectrally linear 
effects (due to aerosols 0 £ other absorbers), but it will be in 
error if the effects are spectrally non-linear (Fig. 13). 



Figure 13. Atmospheric Aerosol and Total Ozone Measurement 


Experimental measurements of atmospheric extinction indicate 
that sometimes there are significant spectral nonlinearities, 
with g 2 ~ 4 X 10"^, and consequent errors in of 5%. There is 

no firm agreement on this problem, since the experimental 
determination of g 2 is not very accurate, and theoretical calcula- 
tions indicate that rather unusual types of aerosols are needed to 
produce the nonlinearities. The effect might well be due to an 
absorbing species. 
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One species which does absorb in the uv is NO 2 . Fortunately, its 
absorption spectrum is reasonably linear so that its error component 
is largely eliminated by the double wavelength pair measurement. 
For a single wavelength pair, though, errors of many percent can arise 
where NO 2 concentrations are high, e.g. 40 ppb averaged over an 
atmospheric column. Typical values are 1-10 ppb, according to 
Basher. Further study is needed to assess the spectral non-linear 
extinction of aerosols and atmospheric trace species, as well as the 
long-term variability of the spectral character. 

Basher then discussed a model of stray light (Fig. 14) that 
indicates a possibility of errors of several percent in the Dobson 
instrument, "always an underestimate." The error arises from the 
interaction of instrument and atmosphere, and is dependent on the 
instrument's sensitivity to the stray light, as well as the solar 
zenith angle and the atmosphere's ozone content. The instrument 
dependence can change slowly with time so a regular monitoring or 
correction of stray light is needed to avoid giving the appearance 
of long-term change in ozone. However, Basher concluded, the 
dependence of the stray light error on zenith angle and ozone itself 
simply adds an extra amplitude to the already existing yearly cycle 
term and is thus of no consequence. 

Examining the case of calibration errors. Basher noted that 
calibration basically consists of ensuring constancy in the wave- 
lengths of the bands and determining each band-pair's relative 
spectral output in the absence of the atmosphere, that is, deter- 
minining extraterrestrial constants. The first part involves simply 
a comparison with some stable wavelength standard. From the point 
of view of trend detection it is desirable to do this at least 
once a year. Quite rough wavelength calibrations would even be 
acceptable if they were made frequently and the "roughness" were 
random. The second part of the calibration is probably the most 
difficult aspect of ground-based total ozone measurement. For such 
measurements the "absence of the atmosphere" can only be gained 
indirectly, and the field measurements involved are time consuming 
and give inconsistent results that reflect as much as several percent 
variation in ozone measurements. 

Basher reinforced Komhyr's point that the solar near-ultraviolet 
relative spectral irradiance may not be constant. For trend analysis, 
extraterrestrial constants need to be constant but not necessarily 
accurate. The instrument component of the constants can be cali- 
brated with standard lamps, although past experience has shown this 
to be no easy task. The measurement of the solar uv component of 
the extraterrestrial constants requires an extra -atmospheric instrument 
location. The stability of this calibrating instrument must be 
better than the suggested 0.3% per year sensitivity of the ozone 
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Figure 14. Model of Stray Light Error in Dobson Total Ozone 
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trend detection method; that is, it must be stable to less than 
0.2% per year in relative spectral response. Moreover, there 
must exist means to independently verify the stability on a 
continual basis. Basher said, "The fact that at present such 
long-term stability is very difficult to achieve, even in the 
laboratory, should not completely discourage the consideration 
of preliminary remote space experiments." 

In conclusion. Basher noted, "Our present lack of knowledge 
of the effects of aerosol attenuation and of other absorbers on 
the ozone measurement and our lack of certainty about the extra- 
terrestrial constants' constancy severely limit our ability to 
interpret measured trends as real ozone variation. It is quite 
possible that aerosols have no significant effect on X /\0 total 
ozone, that the uv solar spectrum is particularly constant, and 
that the 0.3% per year detectability limit is meaningful. However, 
until we can prove these things we must accept a large uncertainty 
in trend detection results, probably more than 1% per year. Of 
course a good deal more study of systematic errors is needed." 

Donald Heath next discussed satellite measurements of ozone 
using the nadir-looking BUV (Backscattered Ultra-Violet) instru- 
ment on NIMBUS 4, which was launched in 1970. This instrument is 
producing data at what Heath termed a cost of "about a dollar 
each over the life of the satellite, including original cost" 

(Fig. 15). NIMBUS G planned for 1978 as well as the TIROS N 
series may carry similar instruments. 

The BUV instrument measures the ultraviolet solar radiation 
which is backscattered by the earth and its atmosphere and compares 
this, at varying wavelengths between 2550 and 3400 angstroms, to 
the incoming uv radiation (Fig. 16). The data are used to estimate 
both vertical profiles of ozone down to the 0.2 millibar level and 
total ozone. The data inversion technique requires a statistical 
quantity for P*, the pressure height of the maximum ozone level as 
a function of latitude and season. Balloon and rocket data help 
provide upper and lower first guesses (Fig. 17). 

So far total ozone data have been released for the first year 
of the satellite's mission and these data are currently undergoing 
a revision based on improved calibration of the instrument. 

"We need guidance to determine long term instrument perfor- 
mance," said Heath. He then went on to discuss the puzzle of know- 
ing which Dobson instruments to use for ground truth. Integrating 
data over the whole globe, he compared BUV and Dobson measurements. 
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Figure 16. Effective Scattering Levels for Solar 
Radiation Scattered in the Nadir Direction of the 
Satellite for all Orders of Scattering (Solar 
zenith angle = 60°; Total ozone = 336 m atm-cm) 
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Figure 17. Balloon and Rocket Data 
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Figures 18 and 19 display two data sets from satellite 
observations. The massive solar proton event of August 4, 1972, 
is seen dramatically in a cavity in the ozone formed over the 
North Pole that persisted for at least two months. This same 
effect was not observed in the southern hemisphere, probably 
because it was winter and the ozone variability was very high. 

Heath concluded that NASA's primary interest is in the 
stable high altitude observations where any effect of CFMs 
should show up. Figure 20 shows NASA's March 1977 assessment 
of the percent change in ozone due to CFMs. 
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Figure 20. Depletion of Ozone due to CFMs, NASA 
Assessment, March 1977 
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John Gille briefly described his Limb Radiance Infrared 
Radiometer (LRIR) results with a special note on the difference 
between accuracy and precision. The LRIR receives radiation 
from a 200- to 300-km-long path through the atmosphere. Cold 
space behind the earth's limb eliminates any problems with back- 
ground radiation and the ray path eliminates interference from 
any lower atmospheric regions. The LRIR produces vertical ozone 
distributions from 15 to 55 kms, as well as temperature distributions 
(15 to 70 km). Water vapor profiles are expected in the future. 

Gille stated that the LRIR accuracy is comparable to that of 
other instruments including the optical rocketsonde, the chemi- 
luminescent sonde, and the balloon ozonesonde. As to precision, 
he noted that the LRIR instrument could get a reading every 25 km 
or 12 to 15 profiles every 4° of latitude. The standard deviation 
for a single profile is a function of altitude, but ranges from 
about 1% at 30 km to 3% at 4.0 km and 10% at 48 km. He noted that 
comparisons between LRIR and Limb Infrared Monitor of the 
Stratosphere (LIMS), to be launched in 1978, would provide a very 
sensitive test of variations at 40 km. For accuracy, he noted 
that the rms agreement with the chemiluminescent rocketsonde 
was 10%, about the stated accuracy of the rocket measurement. 

James Lovill described the Air Force defense meteorological 
polar orbiting satellite which can send 68,400 observations a day 
to the Satellite Ozone Analysis Center (SOAC). Twenty-three 
Dobson stations in 13 countries have agreed to provide special 
ground truth measurements to the SOAC to use with the satellite data. 

Detailing the measurement procedure, Lovill explained the 
instrument, its flow diagram, some specific instrument parameters, 
and the sensor scan geometry (Fig. 21). He then presented satel- 
lite readings (Fig. 22), contoured them over the globe (Fig. 23), 
and provided Arosa readings for July 1977 (Fig. 24) as an example 
of the ground truth measurements. He concluded his presentation 
with a display of correlation coefficients vs. distance for two 
sets of total ozone stations (Figs. 25 and 26). 
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Figure 21. Optical Schematic of Ozone, Temperature and Water Vapor Sensor 
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Figure 24. Arosa Readings for July 1977 
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Figure 25. North American Total Ozone Stations 
(Data Period: 1960-1975) 
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Figure 26. European Total Ozone Stations 
(Data Period: 1960-1975) 


34 




